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We discuss universal scaling properties of (2+l)-flavor QCD in the vicinity of the chiral 
phase transition at vanishing as well as non- vanishing light quark chemical potential (/xz). 
We provide evidence for O(N) scaling of the chiral order parameter in (2+l)-flavor QCD 
and show that the scaling analysis of its derivative with respect to the light quark chemical 
potential provides a unique approach to the determination of the curvature of the chiral 
phase transition line in the vicinity of fii/T — 0. 

§1. Introduction 

It is well understood that the phase structure of QCD at non-zero temperature 
crucially depends on the quark mass values ED'S) While in the case of three degener- 
ate quark masses (3-flavor QCD) a first order transition occurs for sufficiently small 
values of the quark mass, it is expected that the transition in 2-flavor QCD is con- 
tinuous and belongs to the universality class of 0(4) symmetric, 3-dimensional spin 
models. Several studies of the QCD phase diagram as function of two degenerate 
light (m u ,md) and a strange (m s ) quark mass suggest that the region around the 
3-flavor chiral limit, where the QCD transition becomes first order, is indeed small 
and does not include the point of physical light and strange quark masses. Our 

current understanding of the location of 
the region of I s ' order transitions that 
is separated by a line of 2 nd order tran- 
sitions from the crossover region and its 
positioning relative to the physical point 
is presented in the Columbia plot shown 
in Fig. [TJ 

While numerical calculations in 3- 
ffavor QCD gave evidence for the exis- 
tence of a first order transition, many 
of the details of the transition in 2- 
or (2+l)-flavor QCD with light up and 
down quarks are still poorly constrained 
through lattice calculations. Earlier at- 
tempts to verify 0(4) scaling in numer- 
ical studies of 2-flavor QCD using stan- 
dard staggered fermions were not very 
successful in establishing the expected 
universal scaling properties.^ 
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Fig. 1. The phase diagram of 3-flavor QCD. 
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We present here a new analysis of scaling properties using 0(a 2 ) improved stag- 
gered fermion and gauge actions.™ We perform this analysis for (2+l)-fiavor QCD. 
Implicitly we thus also analyze whether the region of first order chiral transitions 
ends in a tri-critical point for a strange quark mass value below its physical value 
or above. We will use the scaling analysis of the chiral order parameter to deter- 
mine the curvature of the chiral phase transition line at non- vanishing values of the 
light quark chemical potential ///, i.e. the first (Taylor) expansion coefficient of the 
transition line in powers of (fn/T) 2 around Hi/T = 0. 

§2. O(N) scaling of the chiral condensate 

In the vicinity of a critical point regular contributions to the logarithm of the 
partition function become negligible. The universal critical behavior of the order 
parameter M of, e.g. 3-dimensional O(N) spin models, is then controlled by a 
scaling function f G that arises from the singular part of the logarithm of the partition 
function, 

M(t,h) = h l ' & f G {z), (2-1) 

with z = t/h}'° and scaling variables t and h that are related to the temperature, 
T, and the symmetry breaking (magnetic) field, H, 

1 T - T c , H , . 

t = T^- 1 , h = — . (2-2) 

to T c ho 

Here j3 and 5 are critical exponents, unique for the universality class of the second 
order phase transition which the system undergoes in the limit (t,h) — > (0, 0). The 
form of the scaling function fc{ z ) is we U known from numerical simulations of 3- 
dimensional 0{N) symmetric spin models.® 

In QCD symmetry breaking arises due to a non-vanishing light quark mass, 
H = rrii/m s , and the corresponding order parameter is the chiral condensate, which 
we write as 

M = M b = NfmsiH'h > (2-3) 

where m s denotes the strange quark mass in lattice units and iV T is the temporal 
extent of the 4-dimensional lattice, N% x N T . One may improve the operator Mb by 
subtracting a fraction of the strange quark condensate.^* This eliminates additive 
linear divergent terms proportional to the light quark masses. Such a subtraction 
is, in fact, mandatory, if one wants to take the continuum limit at finite quark mass 
before taking the chiral limit. This ordering of limits, indeed is needed in order to 
recover the correct 0(4) scaling behavior from calculations performed with staggered 
fermions, which only preserve a global 0{2) symmetry for any non-zero value of the 
lattice spacing. We will be less ambitious here and discuss the chiral limit on lattices 
with fixed temporal extent, N T = 4. In this case, we can only expect to find 0(2) 
rather than 0(4) scaling behavior. However, the scaling functions, fdz), are very 
similar for both universality classes and it thus will be difficult to distinguish 0(2) 
and 0(4) scaling through an analysis of the order parameter alone. Moreover, given 
the large scaling violations observed in earlier studies with staggered fermions,™ 
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already the observation of scaling in terms of a generic 0(N) scaling function at 
non-zero values of the lattice cut-off is a major step forward. 

We show in Fig. [^left) results from a calculation of chiral condensates in (2+1)- 
flavor QCD. The bare strange quark mass (m s ) has been chosen such that the physi- 



cal value of the strange pseudo-scalar mass, m sS = ./2m| — m^, is reproduced. The 
light quark mass has been varied in a range 1/80 < mi/m s < 2/5, which for the light 
pseudo-scalar Goldstone meson corresponds to a regime 75MeV<m ps <420MeV. The 
lattice size has been varied from 16 3 x 4 for the heavier quark masses to 32 3 x 4 for 
the lightest quark masses.™ This insures that finite volume effects remain small in 
the entire light quark mass regime, i.e. in units of the spatial extent N a we always 
have nipsNo- > 3. 

From Fig. [^left) it is obvious that the chiral condensate scales with the square 
root of the quark mass in the low temperature, chiral symmetry broken phase, 



V mo V m s 



This is characteristic for Goldstone- modes in three dimensional 0(N) symmetric spin 
models. In fact, this also is the dominant term characterizing the scaling function 
fc( z ) i n the symmetry broken phase, i.e. for z < 0.™ 

Results for the chiral condensate may be put on the universal scaling curve by 
using the reduced temperature and rescaled symmetry breaking field introduced in 
Eq. 12-21 Of course, this is expected to be possible only close to criticality where 
contributions from regular terms and corrections to scaling are small. The scal- 
ing analysis shown in Fig. fright) has therefore only been performed for the three 
lightest quark mass values, mi/m s < 1/20 and for temperatures close to T c . From 
this one determines the three free parameters, to, ho and T c . As expected results 
for heavier quarks, also shown in Fig. fright), show deviations from the universal 
scaling behavior. Contributions from corrections to scaling become significant for 
mi/m s >l/5, i.e. m ps >300 MeV. This is in contrast to calculations with Wilson 
fermions, where indications for 0{4) scaling have been reported for even large values 
of the pseudo-scalar meson mass.^ 

The scaling analysis of the order parameter provides two non-universal param- 
eters that are unique for QCD, the chiral phase transition temperature, T c , and the 
scale parameter zq = h /to- In the continuum limit both quantities are functions 
of the strange quark mass only. Of course, T c as well as zq are cut-off dependent and 
a proper continuum extrapolation is needed to extract their values in the continuum 
limit. From our analysis on lattices with temporal extent N T = 4 we find zq ~ 7.5. 
A preliminary analysis on lattices with temporal extent N T = 8 suggests that this 
value drops by almost a factor 2.^ A more detailed analysis of the approach to the 
continuum limit thus is needed. We stress, however, that zq is a physical parame- 
ter of QCD. It gives the slope of the quark mass dependence of the pseudo-critical 
temperature, which can be determined from the location of a peak in the chiral 
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Fig. 2. The light quark chiral condensate in lattice units versus the ratio of the square root of the 
light and strange quark masses (left) and its scaling form (right). The right hand figure has 
been obtained using fits for mi/m s < 1/20 and \T — T c \/T c < 1.03 only. 
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/*(*) = - 6 (fa(z) ~ faz) 



(2-5) 
(2-6) 



The scaling functions /q(z) = dfc/dz and f x ( z ) are shown in the right hand part 
of Fig. [H The scaling function f x (z) has a maximum at z p . The dependence of the 
pseudo-critical temperature, T p , on the quark mass is given by the condition that 

z = t/h 1 ^ s = z p , i.e. 

TpjH) - T c = Zp H i//35 _ ^ 2 . 7 ^ 

T c z 

For the 3-d 0(2) universality class the peak in the chiral susceptibility is located at 
z p ~ 1.56. Using this and expressing the symmetry breaking field H in terms of pion 
and kaon masses rather than quark masses, H = mi/m s — 0.52 (m ps /mx) 2 , we find 
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This allows to estimate the phase transition temperature in the chiral limit. With 
m n /m,K ~ 0.27 one finds with our current estimates for zq that the transition tem- 
perture in the chiral limit is about 7% smaller than the crossover temperature at 
physical values of the quark masses. Of course, this still needs to be analyzed closer 
to the continuum limit. 

§3. Curvature of the critical line in the T-fiB phase diagram 



For non- vanishing light quark chemical potential, \i\, the second order chiral 
phase transition persists to exist in the T-[i\ plane. For small values of the chemical 
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potential, /^/T>0, the curvature of the phase transition line can be determined by 
making use of the scaling analysis of the order parameter. In this case the reduced 
temperature variable, t, also depends on the chemical potential as it couples to the 
quark number operator, which does not break chiral symmetry. To leading order it 
contributes quadratically, 
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Fig. 3. The curvature of the transition curve 
in the T-\xb plane. 
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The condition for criticality, t = 0, fixes 
the shape of the transition line for small 
values of m/T. In the scaling regime the 
order parameter M depends on temper- 
ature and quark mass only through the 
scaling variable z. The derivative of M 
with respect to T thus is, up to a con- 
stant, identical to the second derivative 
of M with respect to the chemical po- 
tential hi. This derivative too is related 
to the scaling function f@(z) which we 
have introduced in Eq. 12-61 and which is 
shown in Fig. fright), 
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We note that M<i has properties similar to the chiral susceptibility. It diverges in 
the chiral limit at z = and the peak position in /^(z) can be used to define a 
pseudo-critical temperature at non-zero values of the light quark mass. 

Once the scale parameters to, ho, T c , needed to project the chiral order parame- 
ter onto the 0(N) scaling curve, are known, we can use this information to determine 
the curvature of the critical line, «„, by calculating Mi and by matching the scaling 
curve to the data using an appropriate scaling factor k^. Such a scaling analysis is 
shown in Fig. [Deleft). From this we find for the curvature of the phase transition line 
in the chiral limit the preliminary result k^ ~ 0.035P* which is in good agreement 
with earlier determinations of the curvature of the pseudo-critical line at non-zero 
values of the quark mass performed with different numbers of flavors™ We note that 
the current analysis suggests that the curvature of the transition line, expressed in 
terms oQ/ig ~ 3/U/, is smaller than the experimentally determined freeze-out curve, 
which is well parametrized byp^J 1 
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*' Some caution is needed when translating fii to fis and comparing lattice results with exper- 
imental findings. Lattice calculations have been performed for vanishing strange quark chemical 
potential (in some cases even for 2-flavor QCD), while the freeze-out conditions in a heavy ion 
collision refer to a system with \x s j '/4b — 0.16.^* 
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Fig. 4. The second derivative of the order parameter with respect to the light quark chemical 
potential (left). The 0(2) scaling curve f'c{z), shown together with / x (z) in the right hand part 
of the figure, has been rescaled to match the data. 

§4. Conclusions 

We have discussed universal properties of the chiral condensate of (2+l)-flavor 
QCD in the limit of vanishing light quark masses. It agrees well with expected 0(N) 
scaling predictions. We showed how this can be used to calculate the curvature of 
the second order chiral phase transition line in the vicinity of i^b/T = 0. 

At present this analysis is limited to rather coarse lattices. Calculations closer 
to the continuum limit are needed to get control over the scale parameter zq and the 
phase transition temperature in the chiral limit. 

Acknowledgements 

We thank the organizers of the workshop 'New Frontiers in QCD 2010' at the 
Yukawa Institute for Theoretical Physics, Kyoto, for a very stimulating workshop 
program and for support. This work has been supported in part by contract DE- 
AC02-98CH10886 with the U.S. Department of Energy. 



References 

1) R. Pisarski and F. Wilczek, Phys. Rev. D 29, 338 (1984). 

2) F. R. Brown et al, Phys. Rev. Lett. 65, 2491 (1990). 

3) F. Karsch and E. Laermann, Phys. Rev. D 50, 6954 (1994); 

S. Aoki et al. (JLQCD Collaboration), Phys. Rev. D 57, 3910 (1998); 
C. W. Bernard et al, Phys. Rev. D 61, 054503 (2000). 

4) S. Ejiri et al, Phys. Rev. D 80, 094505 (2009) 

5) J. Engels, S. Holtmann, T. Mendes and T. Schulze, Phys. Lett. B 492, 219 (2000); 
J. Engels and T. Mendes, Nucl. Phys. B 572, 289 (2000). 

6) Y. Iwasaki, K. Kanaya, S. Kaya and T. Yoshi'e, Phys. Rev. Lett. 78, 179 (1997); 
A. Ali Khan et al. (CP-PACS Collaboration), Phys. Rev. D 63, 034502 (2001). 

7) W. Unger (RBC-Bielefeld Collaboration), poster presented at Strong and Electroweak 
Matter 2010, Montreal, Canada, June 29- July 2, 2010. 

8) RBC-Bielefeld Collaboration, in preparation 

9) for a summary of earlier results on the curvature of the transition line see: 
O. Philipsen, Prog. Theor. Phys. Suppl. 174, 206 (2008). 

10) J. Cleymans, H. Oeschler, K. Redlich and S. Wheaton, Phys. Rev. C 73, 034905 (2006). 



